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Abstract: We report results of the performance evaluation of a new hardware correlator in Korea, the Dae¬ 
jeon correlator, developed by the Korea Astronomy and Space Science Institute (KASI) and the National 
Astronomical Observatory of Japan (NAOJ). We conducted Very Long Baseline Interferometry (VLBI) 
observations at 22 GHz with the Korean VLBI Network (KVN) in Korea and the VLBI Exploration of 
Radio Astrometry (VERA) in Japan, and correlated the aquired data with the Daejeon correlator. Eor 
evaluating the performance of the new hardware correlator, we compared the correlation outputs from 
the Daejeon correlator for KVN observations with those from a software correlator, the Distributed EX 
(DiFX). We investigated the correlated flux densities and brightness distributions of extragalactic compact 
radio sources. The comparison of the two correlator outputs show that they are consistent with each other 
within < 8%, which is comparable with the amplitude calibration uncertainties of KVN observations at 
22 GHz. We also found that the 8% difference in flux density is caused mainly by (a) the difference in the 
way of fringe phase tracking between the DiFX software correlator and the Daejeon hardware correlator, 
and (b) an unusual pattern (a double-layer pattern) of the amplitude correlation output from the Daejeon 
correlator. The visibility amplitude loss by the double-layer pattern is as small as 3%. We conclude that 
the new hardware correlator produces reasonable correlation outputs for continuum observations, which 
are consistent with the outputs from the DiFX software correlator. 

Key words: Techniques: interferometric — Instrumentation: interferometers — Radio continuum: galax¬ 
ies — Masers 


1. Introduction 


Very Long Baseline Interferometry (VLBI) is an as¬ 
tronomical observing technique developed in the 1960s 
for measuring the accurate position of compact radio 
sources and obtaining their sky brightness distribution 
at high angular resolution by detecting fringes of noise 
signals arriving at two, or more, radio telescopes from 
the celestial coinpact radio sources (iR.vle fc HewishI 
119601: IClark et. allll967t iMoran et. allll967ll . The noise 
signals arriving at each radio telescope are filtered, 
down-converted, and digitally sampled with being ac¬ 
curately time-tagged by a frequency standard located 
at each observatory. The sampled signals are recorded 
using magnetic tape sy stems such as the Mark I sys- 
tem ( Bare et. all llOOlt . the Mark H system (GlarkI 
Il973ll . the Mark HI system ( Rogers et. all 119^ 1. the 
Mark IV system (IWhitnevI Il993ll . and the S2 sys¬ 
tem ( Wietfeldt et. all 199611 . or to hard d isk systems 
such as the Mark 5 sy stem (|Whitnevll2002ll and the K5 
system (iKondo et. all 1200^ . The recording media are 
shipped to a correlator center and played back. At that 
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moment, one of the main roles of the correlator is to 
detect a fringe between the two signals from radio tele¬ 
scopes after correcting for geometric and instrumental 
delays of the signals. Pre-estimates of the delays, apri- 
ori data are applied to the correlation of the signals 
from the two radio telescopes (or from a baseline). The 
correlation output in this stage is known as the visi¬ 
bility. The correlation output streams are formatted 
in Flexible Image Transport System ( FITS) interfer - 
ometry Data Interchange Gonvention (jGreise 3|2Q0i). 
Through a post-correlation process, the visibilities are 
used for measuring the accurate positions of celestial 
radio sources and recovering their sky brightness dis¬ 
tributions. The post-correlation process is performed 
generally with astronomical processing programs: for 
example, the Astronomical Image Processing System 
(AIPS). 


A number of VLBI correlators have been devel¬ 
oped to detect the fringes from VLBI observations. 
The VLBI correlators include hardware and soft- 


ware correlators (iBare et. all 119671: Moran et. all 

1967 

Naoier et. all 11994: Wilson et. al 119961; Casse 

1999 

Gar Ison 

19991; Horiuchi et al. 

20001: Deller et. all 

2007 
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[mil). Technical progress in parallel computing and 
high-speed networks enables the construction of soft¬ 
ware correlators at low price and with short develop¬ 
ment perio d. Especially software c orrelator such as 
the DiFX (iDeller et. all l2007l l201ll) became popular 
because of its easy installation and extensive support 
from its user groups. The price per baseline and price 
per unit data bandwidth are almost the same in soft¬ 
ware correlators, while the prices decrease significantly 
with the number of baseline and data rate in hardware 
correlators. Therefore, hardware correlators become at¬ 
tractive when the target VLBI array consists of a large 
number of antennas with a wide bandwith system. In 
addition, the electric power consumption of hardware 
correlators is much less than that of software corre¬ 
lators, significantly reducing the total operation bud¬ 
get. These are the main reasons why connected arrays, 
which usually have more elements and wider bandwidth 
than VLBI arrays, still prefer hardware correlators. 

A new hardware correlator (the Daejeon correlator, 
FigurejT]) was developed in 2006-2009 by the Korea As¬ 
tronomy and Space Science Institute (KASI) and the 
National Astronomical Observatory of Japan (NAOJ). 
The Daejeon correlator will be the main correlator for 
the East Asian VLBI Network (EAVN) co nsisting of 
the Korean VLBI Network (KVN) in Korea (|Lee et. all 
1201 1L 12014 . the Japanese VLBI Network (JVN) in¬ 
cluding the VLBI Exploration of Radio Astrometry 
(VERA) in Japan, and the Chinese VLBI Network 
(CVN). This is one of the main motivations for develop¬ 
ing the new hardware correlat or rather than adopting a 
software correlator like DiFX ([Deller et. alll200^l201ll1 . 
The capability of efficient correlation for large number 
of stations (e.g., >20 for EAVN stations at 22 GHz) is 
one of the advantages of the hardware correlator. In 
2010 we started commissioning operations of the corre¬ 
lator, and began to evaluate its performance using test 
observations with the combined network of KVN and 
VERA (Figure O. 

The Daejeon hardware correlator is located in the 
Korea-Japan Correlation Center (KJCC), Korea. The 
Daejeon correlator consists of several VLBI data play¬ 
back systems, a Raw VLBI Data Buffer (RVDB), a 
VLBI Correlation Subsystem (VCS), and a data archive 
system. The Daejeon correlator was aimed to corre¬ 
late data obtained from various VLBI networks in East 
Asia: the KVN, the KVN and VERA combined net¬ 
work, and the East Asia VLBI Network. As these 
use different recording systems, the Daejeon correla¬ 
tor has several VLBI data playback systems: Mark 5B, 
VERA2000, OCTADISK, etc. They have different in¬ 
terface for data transmission. The difference in the for¬ 
mats of the playback system led to the introducion of 
the RVDB system, a big data server with several inter¬ 
faces. Data from the RVDB are transferred to the VCS, 
the main part of the Daejeon correlator, and correlated 
with proper control parameters provided by the cor¬ 
relator control and operation computers. The VCS is 
able to process a maximum of 16 stations with a max¬ 
imum recording rate of 8192 Mbps, and 8192 output 



channels for VLBI observation data. The correlation 
results from the VCS are then transferred to the data 
archive system. Detailed description of the correlator 
will be presented elsewhere (C. S. Oh et al. 2015, in 
preparation), and its current status can be found in the 
Internet 13 

Performance of the new correlator can be evaluated 
via investigating the outputs from the whole cycle of 
the VLBI observations: data acquisition, correlation, 
and post-correlation process. Careful comparison of the 
output from the new correlator with that from an ex- 
isting, reliable correlato r (e.g., DiFX as described in 
iDeller et. alll2007l l201l[l should be conducted. 

We conducted test observations of compact radio 
sources using KVN and VERA, correlated the observed 
data with the Daejeon correlator and the DiEX soft¬ 
ware correlator, and performed careful comparison of 
the outputs. An extensive evaluation using the Daejeon 
correlator, the DiFX correlator, and Mitaka FX corre¬ 
lator in NAOJ (IShibata et. all[l998l1 for the KVN and 
VERA observations has been conducted. Eull results 
of the extensive comparison will be reported elsewhere 
(C. S. Oh et al. 2015, in preparation). In this paper, we 
report the results of the comparison between the Dae¬ 
jeon correlator and the DiFX correlator using the ob¬ 
servations only with the KVN. In Section 2, we describe 
the observations, correlation, and the post-correlation 
process. The results of this comparison are reported 
in Section 3, and we discuss an unusual pattern of the 
correlation output in Section 4. We make conclusions 
in Section 5. 

2. Observations and Data Reduction 

2.1. Observations Using KVN and VERA 

We observed compact radio sources, 3C 454.3, 3C 345, 
NRAO 512 and Sagittarius A* (Sgr A*) at 22 GHz, 
on 2011 January 28 UT 19:00 and January 29 UT 
04:00, using KVN and VERA (7 stations): KVN Yon- 

^http://kjcc.kasi.re.kr 
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VERA 


http;//veraserver.mtk.nao.acJp/index.htn 


Antenna Diameter 20m (250 fl m) 
Observing band 2,8,22.43GHz 
Maximum baseline 2273km 
Minimum baseline 1000km 
2 beam phase referencing 
-> Inst. Path error <0.1 mm 


Figure 2. Korean VLBI Network (KVN) operated by KASI, Korea (left) and VLBI Exploration of Radio Astrometry (VERA) 
operated by NAOJ, Japan (right). 


sei (KY), KVN Ulsan (KU), KVN Tamna (KT), Mizu- 
sawa (MIZ), Ishigakigima (ISH), Ogasawara (OGA), 
and Iriki (IRK). In order to evaluate the various cor¬ 
relation performances of the Daejeon correlator, the 
several observing modes and scan lengths were consid¬ 
ered. For this comparison, we used longer scans (10-20 
min) for all sources. The observing frequency is 22.00- 
22.45 GHz with a selected bandwidth of 256 MHz and 
in left circular polarization (LGP). The received sig¬ 
nals within the frequency bandwidth are 2-bit quan¬ 
tized by AD converters located in the telescope cabin 
and transferred to the observing building via optical fi¬ 
bre. The transferred digital signals are divided into 16 
sub-bands (IFs) by a digital filter bank, and recorded 
in magnetic disks or tapes by two data acquisition sys¬ 
tems: the Mark 5B system (disk-based) in KVN and the 
VERA2000 system (tape-based) in VERA. The record¬ 
ing rate is 1024 Mbps. 

2.2. Correlation Using Daejeon Correlator and DiFX 

The recording media were sent to the KJCC. The 
shipped media were played back using Mark 5B for 
KVN data and VERA2000 for VERA data. The data 
from the playback systems were transferred to the 
RVDB to be correlated in the VCS. Finally, the data 
were correlated using a 1.6384 s correlator integration 
time and 128 spectral channels across each 16 MHz 
band. The correlation was performed in 2013 July. 
The data recorded at KVN stations were also corre¬ 
lated with the DiFX software correlator. Each scan 
recorded in the Mark 5B disk packs was extracted by 
using the utility FUSEQ and was transferred to the 
hie server through the 10 Gigabit network. Scans for 
strong sources were Hrst used for searching the clock 
offset and rate of each station, and then all scans were 

^ http://fuse.sourcef orge.net 



Figure 3. Schematic diagram of VLBI data reduction with 
AIPS 


correlated. The correlator integration time was 2.048 s 
for DiFX and the number of spectral channels was 128 
across each 16 MHz band. 

2.3. Post-Correlation Process Using AIPS and DIFMAP 

We performed further processing using AIPS. We fol¬ 
lowed a standard processing procedure as described in 
Figure |31 The correlated output in FITS format was 
uploaded into an AIPS catalog by the FITLD task, 
and information of the FITS hie were indexed with 
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the INDXR task. In using FITLD, we did not apply 
digital correction (Digicor = -1), since it is not nec¬ 
essary for the two correlators. We sorted the indexed 
catalog in an order of time and baseline (TB) and if 
necessary, indexed them again. FITS files from the 
correlator that were divided into several files in time 
were combined using the DBCON task, while those di¬ 
vided in frequency were combined with the VBGLU 
task. At KVN and VERA, we sampled the received 
signals using the digital samplers. During the digital 
sampling, there could be some amplitude errors of cross¬ 
correlation power spectrum. These errors can be cor¬ 
rected using the ACCOR task based on the amplitude 
of the auto-correlation power spectrum for each station. 
These corrections for the amplitudes are stored in a so¬ 
lution (SN) table, separate from the FITS file itself. In 
order to apply the correction in the SN table, we used 
the CLCAL task and generated a new calibration (CL) 
table. At this stage, one may average the FITS data in 
frequency using the AVSPC task. We did not average 
the FITS data for this work. The residuals of the fringe 
delay and delay rate after the correlation were searched 
using the FRINC task. This is known as fringe-fitting. 
FRINC provides us with antenna-based solutions of the 
delay and delay rate based on baseline-based solutions. 
The solution interval for the fringe-fitting should be 
pre-determined according to the coherence time of the 
observation at the observing frequency. For this obser¬ 
vation at 22 CHz, we used 1 min as the solution interval 
of the fringe-fitting. Once we obtained the solutions of 
the delay and delay rate, we calibrated the phase using 
CLCAL. After the phase calibration, we performed the 
amplitude calibration in order to correct for the atmo¬ 
spheric opacity change and for the amplitude errors due 
to atmospheric fluctuations. We used system temper¬ 
atures and antenna gain measured at each observatory 
for converting the correlation coefficient to sky bright¬ 
ness and correcting for the amplitude errors. We stored 
the calibration information in the TY and CC tables us¬ 
ing the ANTAB task and produced the solution of the 
amplitude calibrations using the APCAL task. We did 
not correct for the effect of the bandpass filter on the 
spectrum shape in order to investigate the correlation 
effect on the spectrum shape. We divided the FITS 
data for each source with applying all calibration in¬ 
formation and exported into uv-Sle using FITTP. For 
the purpose of the comparison of images, we used the 
central 115 channels of each IF, excluding the first 6 
and last 7 channels. 

After the phase and amplitude calibration, we made 
the contour maps o f the target sources usi ng the Caltech 
DIFMAP software (|Shepherd et al.l[T994f) . We averaged 
the uv-data, (or visibility data) in a time interval of 30 s, 
and flagged outliers in amplitude in order to increase 
the data processing speed and decrease random error of 
individual visibility data. We should note that the aver¬ 
aging time interval of 30 s is shorter than the coherence 
time of this observation. As a first step of mapping (or 
imaging) with DIFMAP, we fitted a point-source model 
to the visibility and self-calibrated the phase according 


to the model, in order to find the converged model to 
the visibility. As a second step, we used the CLEAN 
deconvolution algorithm and the amplitude- and phase- 
self-calibration alternatively. The CLEAN deconvolu¬ 
tion is a technique to deconvolve the visibility data to 
find the true sky brightness distribution of the target 
source by establishing a group of delta-function models. 
The self-calibration is an algorithm to reduce the differ¬ 
ence between the models and visibility phase/amplitude 
using closure phase and closure amplitude. We should 
note that the amplitude-self-calibration can not be per¬ 
formed for the observations with three or less stations. 
Therefore, for the analysis of these KVN data, we did 
not perform amplitude-self-calibration. We evaluated 
the quality of the final map, by investig ating the resid¬ 
ual no ise in the image as described in iLobanov et. all 
HOO^. We can quantitatively estimate the noise in the 
final CLEANed image based on the ratio of image noise 
rms to its mathematical expectation, ^r- Suppose that 
a residual image has an rms Ur and a maximum abso¬ 
lute flux density |sr|. For Caussian noise with a zero 
mean, the expectation of Sr is 


= CTr 


72 In ( 


Nr, 


V / 


1/2 


( 1 ) 


where Vpix is the total number of pixels in the image. 
The ratio is given by 


7 — ■Sr/Sr,exp- (2) 

When the residual noise is similar to Caussian noise, 
t I. If > 1, not all the structure has been ade¬ 
quately cleaned or recovered; if ^r < 1, the image model 
obtained has an excessively large number of degrees 
of freedom. In order to compare quantitative infor¬ 
mation of the final images, we used circular Caussian- 
component models to fit the self-calibrated data, yield¬ 
ing the following modelfit parameters: the total flux 
density ^t, size d, radial distance r, and position angle 
9 of each component. 


3. Results 

3.1. Visibility Comparison with DiFX 

After the calibration and data reduction, we compared 
the calibrated output of the Daejeon correlator with 
that of DiFX. For a careful comparison, we used the 
same parameters and data reduction procedures for the 
two cases, except for the correlator integration time 
(1.6384 s for Daejeon and 2.048 s for DiFX) and the 
delay models (a Mitaka delay model for Daejeon corre¬ 
lator and a CALC9 delay model for DiFX). 

Figures 0115] show the visibility amplitudes on the 
sources, 3C 454.3, 3C 345, and NRAO 512 for IF 2- 
14 and three KVN baselines from both correlators as 
a function of time. The whole frequency channel (128) 
data have been vector-averaged. We fitted the visibility 
amplitudes for each baseline with a first-order polyno¬ 
mial V{t) = At+ B using a linear least-squares method. 
The fitted models and averaged values are shown in 
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Figure 4. Comparison of the visibility amplitude as a function of time for the source 3C 454.3 as calculated by the DiFX (red 
dots) and the Daejeon correlators (blue dots) for IF 2-14 and three KVN baselines: KU-KT, KY-KU, and KT-KY during 
9 min after 2011 January 29 04:20:00 UT. The visibility amplitude is in Jy and the time is in seconds from the start of the 
scan. The solid lines represent linear fits to the visibility amplitude for the DiFX (red line) and the Deajeon correlator (blue 
line). The numbers on each panel indicate the averaged amplitude for each correlator over the time period. 


Figures miHl The averaged visibility amplitudes over 
IF 2-14 within full frequency channels and central 115 
channels for individual baselines and sources are sum¬ 
marized in Table [TJ In Table [TJ we summarized the 
mean visibility amplitudes of IF 2-14 for two cases of 
using full frequency channels in IF and of using only 
central 115 channels. The visibility amplitudes of the 
target sources are in the range of 0.5-22 Jy. The vis¬ 
ibility amplitudes for the Daejeon correlator seem rel¬ 
atively lower than those for DiFX. We found that the 
averaged visibility amplitudes of individual baselines for 
the Daejeon correlator are lower by < 8% than those of 
DiFX for all sources and all baselines. 

The amplitude difference between two correlators can 
be attributed mainly to an unusual pattern of the vis¬ 
ibility amplitude of the Daejeon correlator, so called 


double-layer pattern, as clearly shown in Figured The 
double-layer pattern consists of two main patterns at 
higher and lower amplitude values. One pattern is sim¬ 
ilar to that of DiFX (normal pattern). The other has 
lower values in amplitude than the main pattern and 
sometimes is similar to a sinusoidal pattern. For the 
baseline KU-KT in IF 2 in Figure 01 one third of the 
visibility data are along the level of about 16 Jy, and 
the remaining data are at about 18 Jy. For the baseline 
KU-KT in IF 7, one third of the visibility data show 
a sinusoidal pattern with its amplitude of about 2 Jy, 
that is about 15% of the flux density. The double-layer 
patterns appear different in IF (i.e., in frequency) and 
are present in all baselines and IF bands except for IF 
8 and 9. Among the IF 2-14 for the baseline KU-KT, 
for IF 2-5 and IF 12-15 the double-layer patterns are 
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Figure 5. Same as Figure[4]but for the source 3C 345 and the time range during 19 min after 2011 January 28 19:30:00 UT. 


shown as two parallel patterns, and for IF 6-7 and IF 
10-11 the double-layer patterns appear like the sinu¬ 
soidal pattern. The visibility amplitudes for IF 8 and 9 
do not show such a double-layer pattern. It seems that 
the period of the sinusoidal pattern varies in time and 
frequency. The periods for IF 6 and 11 seems twice 
smaller than those for IF 7 and 10, respectively. We 
found also that there is a rough symmetric trend of 
the double-layer patterns between IF 2-8 and IF 9-15. 
This implies that the double-layer patterns may hap¬ 
pen in a correlation stage related with frequency. The 
double-layer patterns are clearer for the observations 
of stronger sources like 3C 454.3 with a correlated flux 
density of ~20 Jy (FigureHj), whereas they are less clear 
for relatively weaker sources like 3C 345 and NRAO 512 
with flux densities of 1-4 Jy (Figures [S] and H]). This 
is because the amplitude of the double-layer pattern is 
about 15% and it is similar to the amplitude calibra- 
tion uncertainty o f KVN observati ons at 22 GHz (see 
iPetrov et. alll20l3 iLee et. a]ll2ni4ll . However, the 8%- 


difference of visibility amplitudes between DiFX and 
the Daejeon correlator appears in all sources. 

Figure [7] shows the visibility phases of the source 
3C 454.3 for IF 2-14 and three KVN baselines from 
both correlators as a function of time. The whole fre¬ 
quency channel (128) data have been vector-averaged. 
We found that the visibility phases show small differ¬ 
ences due to the different delay models used for the two 
correlators. 

Figure [5] shows a comparison of the cross-correlated 
spectrum, i.e., the visibility amplitude and phases as a 
function of frequency in one 16 MHz band, IF 2. The 
spectra have been vector-averaged over one minute time 
range. The phases of the cross-correlated spectrum are 
consistent with each other within <2 degrees, whereas 
the amplitudes are different by about 10% on average. 
The amplitude difference varies in frequency. At the 
band edges, the difference gets larger, whereas at the 
band center, it becomes a bit smaller. This indicates 
that the spectrum shape of Daejeon correlator output 
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Figure 6. Same as Figure[4]but for the source NRAO 512 and the time range during 19 min after 2011 January 28 19:10:00 
UT. 


Table 1 

Comparison of visibility amplitudes 


Source 

(1) 

Correlator 

(2) 

Sku- 

-KT 

Sky- 

-KU 

Skt- 

-KY 

FULL BW 
(3) 

115CH 

(4) 

FULL BW 
(5) 

115CH 

(6) 

FULL BW 
(7) 

115CH 

(8) 

3C 454.3 

DiFX 

20.2 ± 0.547 

21.2 ±0.576 

21.5 ±0.392 

22.6 ±0.413 

21.0 ±0.554 

22.0 ±0.586 


Daejeon 

18.8 ±1.04 

19.8 ±1.09 

20.0 ±1.07 

21.1 ± 1.13 

19.4 ±1.17 

20.4 ± 1.23 


5'Daejeon/“SDiFX 

0.93(1.9) 

0.93(1.9) 

0.93(2.7) 

0.93(2.7) 

0.92(2.1) 

0.93(2.1) 

3C 345 

DiFX 

3.21 ±0.169 

3.37 ±0.181 

4.37 ±0.158 

4.59 ±0.170 

3.26 ±0.167 

3.43 ±0.179 


Daejeon 

2.93 ±0.236 

3.09 ±0.252 

3.99 ±0.291 

4.21 ± 0.308 

2.97 ±0.239 

3.13 ±0.255 


‘S'Daejeon/iSDiFX 

0.91(1.3) 

0.92(1.4) 

0.91(1.84) 

0.92(1.8) 

0.91(1.4) 

0.91(1.4) 

NRAO 512 

DiFX 

0.861 ±0.163 

0.908 ±0.176 

0.858 ±0.140 

0.903 ±0.153 

0.837 ±0.149 

0.881 ±0.161 


Daejeon 

0.793 ±0.188 

0.839 ± 0.203 

0.792 ±0.169 

0.837 ±0.183 

0.768 ±0.172 

0.812 ±0.186 


“S'Daejeon/‘S'oiFX 

0.92(1.2) 

0.92(1.2) 

0.92(1.2) 

0.93(1.2) 

0.92(1.2) 

0.92(1.2) 


Column designation: 1 - source name; 2 - correlator; 3-9 - mean of visibility amplitude for KU-KT, KY-KU, and KT-KY in Jy, or the 
ratio of the mean values for Daejeon to those for DiFX with the ratio of their standard deviations in parentheses: Columns 3,5,7 for 
values using data in full bandwidth and Columns 4,6,8 for using data excluding first 6 channels and last 7 channels in each band (hence 
averaged over 115 channels). 
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Figure 7. Comparison of visibility phase as a function of time for 3C 454.3 from DiFX (red dots) and the Daejeon correlator 
(blue dots) outputs for IF 2-14 and three KVN baselines: KY-KT, KY-KU, and KT-KY during 9 min after 2011 January 
29 04:20:00 UT. The visibility phase is in degrees and the time is in seconds from the start of the scan. 


is a bit different from that of DiFX. 

3.2. Image Comparison with DiFX 

In addition to the comparison of the visibility ampli¬ 
tudes and phases in time and frequency, we compared 
the imaging results from the two correlator outputs for 
all target sources, 3C 454.3, 3C 345, NRAO 512, and 
SgrA*. As mentioned previously, we used the same 
procedure of imaging for the comparison: (a) fitting 
a point source to the visibility, (b) CLEAN and phase- 
self-calibration, and (c) modelfit with circular Gaussian 
components. In Figure IHl we compare CLEANed im¬ 
ages for the DiFX and Daejeon correlators. We show 
the contour maps of each source for DiFX in the left 
panel, and those for Daejeon correlator in the right 
panel. There are circular Gaussian models on top of 
the contour maps. The x and y axes are in units of 
mas. The sources are 3C 454.3, 3C 345, NRAO 512, 


and SgrA* from top to bottom. For each contour map, 
the source name and the observation date are given in 
the upper left corner of the map. In the lower right 
corner of the map, we show the peak flux density and 
the lowest contour level in units of mJy. The shaded 
ellipse in the lower left corner of the map represents the 
FWHM of the restoring beam of the image. In all of 
the images, the contours are drawn at 1, 1.4, ..., 1.4" 
(logarithmic spacing) of the lowest flux density level. 
In the right panel, we showed the visibility amplitudes 
as a function of the uv radius. The x axis of the plot 
represents the uv radius in units of 10®A, where A is the 
observing wavelength, which is the length of the pro¬ 
jected baseline used to obtain the visibility data. The 
y axis of the plot shows the amplitude of each visibility 
point (i.e., correlated flux density) in units of Jy. The 
quality of the final images was investigated by estimat¬ 
ing the quality of the residual noise, as described in 















































A New Hardware Correlator in Korea: Performance Evaluation with KVN 


9 


10 

a 5 

oc 

3 0 

I -5 

CL 

-10 




.•x'. • '•-'• 


KU-KTIF2 DiFX 

Daejeon 

20 40 60 80 

channels 


*-.*^**.V~.*.*H.*“**V 


100 120 


B 0-9 - 

E 
< 


>• .•'••**^*»V**»»*^»V*»****N,*»,**»N»*». 




Daejeon/DiFX 


60 80 
channels 


10 

a 5 

oc 

3 0 

J -5 

CL 

-10 



•• 


KY-KUIF2 DiFX 

Daejeon 

60 80 
channels 



60 80 
channels 


10 

5 

0 

-5 

-10 



a 10 

00 

3 5 

1 ° 

’■a -5 

-M*,S%«*,.Hhh**^«*««*^**.»***^****^*****«**“^*^»**^**»*^*»*^**»**^******^* “•****• • •*** •- 

** * • , V^*,v.^*.v*«**»«***«x********^*«********•••s***»***»*« ••••, • 

,• s 

£ 

o 

1 

Daejeon/biFX • 

*••••• • • 

• ••• V *“ * ■ • 

• ••• * • ***•**•*•*%•*,• •• 

.. . 

KT-KY IF2 DiFX 

Daejeon 

< 

0.8 

• • ••• /, •* 


60 80 
channels 


60 

channels 


Figure 8. Left: Comparison of the visibility phase (upper panel of each spectrum) and amplitude (lower panel of each 
spectrum) as a function of frequency for 3C 454.3 from the DiFX (red dots) and the Daejeon correlator (blue dots) outputs 
for IF 2: 22.034-22.050 GHz band and three KVN baselines: KU-KT (top), KY-KU (middle), and KT-KY (bottom) during 
1 min after 2011 January 29 04:20:00 UT. The visibility phase is in degrees spanning a range of ±10° and the visibility 
amplitude is in Jy. Right: Phase difference (upper) and amplitude ratio (lower) for two correlators. The phase is shown 
with the same scale as the left panels while the amplitude ratio spans the range of 0.8-1.0. 


Section 12.31 More detailed parameters of the images 
presented in Figure IHl are summarized in Table [2] For 
each image. Table [2] lists the source name, the param¬ 
eters of the restoring beam (the size of the major axis, 


Bmajj the minor axis, and the position angle of 

the beam, Spa), the total flux density, S't, the peak flux 
density, Sp, the off-source rms, u, the dynamic range of 
the image, D, and the quality of the residual noise, 
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for each image. The image results show that the total 
flux and peak flux densities for all sources are consis¬ 
tent within 8% for the two correlator outputs. The 
difference is similar to that of the visibility amplitude 
averaged for each baseline as summarized in Table [T] 
In order to investigate the flux differene in more detail, 
we made images using data for IF 8 and 9, for which 
the double-layer pattern does not appear in Figure |H 
We compared the total and peak flux density for the 
images free of the double-layer-pattern, and found that 
their difference is on average about 5% for all sources, 
as summarized in Columns (11) and (12) of Table [2] 
This may indicate that a part of the total difference 
of 8% is caused by other differences between the two 
correlators that are not related with the double-layer 
pattern. We also compared the flux density, size, and 
position of the core and jet components of all sources 
using circular Gaussian model-fitting. Table [3] lists the 
parameters of each model-fit component: the total flux, 
St, size, d, angular distance from the central compo¬ 
nent, r, and position angle, 6 (the location of the jet 
component with respect to the core component). The 
modelfit parameters show that the total flux densities 
of core and jet components of almost all sources are 
again consistent within 8% for the two correlator out¬ 
puts, whereas the sizes of the Gaussian components are 
different by < 40%. For 3G 454.3 and 3G 345, we were 
able to modelfit the jet components for the two cor¬ 
relator outputs with very consistent positions to each 
other. 

4. Discussion 
4.1. Double-Layer Patterns 

The most prominent, apparent difference between cor¬ 
relation outputs from the Daejeon and DiFX correlators 
is the double-layer patterns as shown in Figure E) The 
double-layer patterns vary in IF (i.e., in frequency) and 
are present in almost all baselines including KVN (and 
also VERA stations). When we investigated the visibil¬ 
ity amplitude of 3G 454.3 for all baselines of KVN and 
VERA, we found that the amplitudes for all baselines 
show the double-layer patterns, except for baselines 
with larger scatter (e.g., Ishigaki-KVN). For VERA 
baselines, the double-layer patterns consist of two par¬ 
allel patterns, and for KVN baselines we see clearly 
sinusoidal patterns, as discussed above. This is consis¬ 
tent with the results for the KVN baselines of fainter 
sources. 

It is important to investigate possible reasons caus¬ 
ing the double-layer pattern and to improve the perfor¬ 
mance of the Daejeon correlator. A detailed investiga¬ 
tion of the reasons behind the double-layer pattern is 
underway by the KJGG engineering team and will be 
reported elsewhere. Here we offer a brief discussion of 
possible ways this pattern can be produced. Possible 
reasons include the following: 

• Fringe rotation error. The plot of the visibil¬ 
ity as a function of time shows very periodic pat¬ 
terns, which may indicate some fundamental prob¬ 


lem in the correlator: for example, the delay shifter 
or the phase rotator is not updated as it should. 
However, as the pattern appears only in the ampli¬ 
tude, it may not be the case. Moreover, same pat¬ 
terns are seen on two separate channels: IF3=IF14, 
IF4=IF13, IF5=IF12, IF6=IF11,IF7=IF10. IF8 
and IF9 are both OK, IF2 has no match, and IFl 
is not plotted. Since these matching channels have 
different frequencies, the fringe rotator is not the 
problem. 

• Timing problem between normalization and 
accumulation. By looking at the plot more 
closely, we found that some integration intervals 
may have lost part of the data. The integration 
time of 1.6384 s for the Daejeon correlator is so 
unusual that may modulate with other time in the 
correlator system. Integer number of seconds of 
integration time, or shorter intervals, would give 
better results. It is possible that the IF pattern 
mentioned above can be compared with the order 
in which the control computer gets the data from 
the hardware. We found that this problem does not 
affect the correlation ouput of DiFX since there is 
no data loss on the normalization and accumula¬ 
tion in DiFX with the time interval of 2.048 s. 

• Problem with the geometry applied in the 

correlator. The period of the sinusoidal pat¬ 
tern of the double-layer patterns seems to vary 
in time and in frequency. This could indicate 
that there is a problem with the geometry ap¬ 
plied in the correlator. The correlation output of 
VLB! observations with a finite bandwidth Av is 
r oc cos{2'KVt)Tg — (j>), where Tg is the geo- 

met ric delay and vg is the observing frequency (see 
e.g.. iThomoson |[i999ll . The amplitude of the cor¬ 
relation output is modulated by a sinc-function en¬ 
velope as . Usually the modulation is well 

modeled by the geometric model used during the 
correlation. However, if there is an error in geom¬ 
etry, i.e.. At = Tg — Tm, where Tm is the geometric 
time delay of the geometric model, then we may 
expect a periodic pattern in the amplitude of the 
correlation output equal to 

case, the sinusoidal pattern may correlate with the 
baseline length. However it is hard to tell wether 
we see the correlation in the KVN observations, 
since the baseline length of KVN is in the range 
of 305-477 km. Moreover, the geometric problem 
should affect the phase of the correlation output 
as cos{2TTVQTg — (f). Since we see no prominent 
difference in phases between Daejeon and DiFX 
correlators, we can exclude this possibility. 

Although the amplitude of the double-layer patterns 
is as large as 10-15% the final effect to the flux density 
in the GLEANed image is less than 8% (maybe 3%, see 
below) because the visibility data are averaged in time 
and frequency for the final imaging. However this effect 
should be investigated more carefully with further test 
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Figure 9. Left: CLEANed images using the DiFX correlator outputs. Circular Gaussian models are on top of the contour 
maps. The axes of each map are the relative R.A. and decl. offsets from the tracking center in milliarcseconds. The lowest 
contour level is shown in the lower right corner of each map. The contours have a logarithmic spacing and are drawn at 1, 
1.4, ..., 1.4" of the lowest contour level. Middle: CLEANed images using the Daejeon correlator outputs in the same format 
as the left panel. Right: distributions of the visibility amplitude for DiFX (red) and the Daejeon correlator (blue) against 
the uv radius. The x axis shows the uv distance in 10®A, and the y axis represents the visibility amplitude (correlated flux 
density) in Jy, averaged over 30 s. Image parameters of each image are summarized in Tabled 
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Table 2 

Comparison of imaging parameters 


Source 

(1) 

Correlator 

(2) 

-^maj 

(3) 

-^min 

(4) 

SpA 

(5) 

5t 

(6) 

(7) 

cr 

(8) 

D 

(9) 

G 

(10) 

S't.iFs-g 

(11) 

•SpaFs-g 

(12) 

3C 454.3 

DiFX 

5.829 

2.937 

-65.1 

24.10 

22.27 

61.53 

362 

0.64 

24.17 

22.28 


Daejeon 

5.835 

2.939 

-65.1 

22.26 

20.71 

57.12 

363 

0.76 

23.23 

21.41 


Daejeon/DiFX 

- 

- 

- 

0.92 

0.93 

0.93 

1.0 

- 

0.96 

0.96 

3C 345 

DiFX 

7.323 

3.111 

-55.1 

4.703 

3.978 

12.14 

328 

0.57 

4.656 

4.018 


Daejeon 

7.333 

3.114 

-55.1 

4.321 

3.640 

15.17 

240 

0.56 

4.390 

3.790 


Daejeon/DiFX 

- 

- 

- 

0.92 

0.92 

1.2 

0.73 

- 

0.94 

0.94 

NRAO 512 

DiFX 

7.741 

3.151 

-54.1 

0.887 

0.871 

2.387 

365 

0.53 

0.897 

0.876 


Daejeon 

7.746 

3.153 

-54.1 

0.818 

0.797 

2.884 

276 

0.53 

0.855 

0.829 


Daejeon/DiFX 

- 

- 

- 

0.92 

0.92 

1.2 

0.76 

- 

0.95 

0.95 

SgrA* 

DiFX 

7.033 

4.493 

-33.5 

1.011 

0.822 

5.023 

164 

0.88 

1.019 

0.828 


Daejeon 

7.042 

4.494 

-33.6 

0.946 

0.772 

4.715 

164 

0.80 

0.988 

0.801 


Daejeon/DiFX 

- 

- 

- 

0.94 

0.94 

0.94 

1.0 

- 

0.97 

0.97 


Column designation: 1 - source name; 2 - correlator; 3 - major axis[mas]; 4 - minor axis[mas]; 5 - position angle of the major axis [°] 
of the restoring beam; 6 - total flux density [Jy]; 7 - peak flux density [Jy beam“'^], and the ratio of the total and peak flux densities 
of Daejeon correlator to those of DiFX; 8 - off-source RMS in the image [Jy beam“^]; 9 - Dynamic range of the image (D = Sp/tr); 
10 - quality of the residual noise in the image; 11-12 - total flux density [Jy] and peak flux density [Jy beam“^] for images using data 
for IF 8 and 9. 

Table 3 

Comparison of model fit parameters 


Source 

(1) 


St 

aejeon 


d 

^Daejeon 


r 


e 

DiFX 

(2) 

Daejeon 

(3) 

DiFX 

(5) 

Daejeon 

(6) 

DiFX 

(8) 

Daejeon 

(9) 

DiFX 

(10) 

Daejeon 

(11) 

'S'dipx 

(4) 

‘^DiFX 

(7) 

3C 454.3 

22.4 

20.8 

0.93 

0.315 

0.199 

0.63 

- 

- 

- 

- 


0.872 

0.835 

0.96 

1.59 

1.48 

0.93 

4.05 

3.83 

-113 

-127 


0.780 

0.624 

0.80 

3.11 

2.73 

0.88 

7.21 

6.14 

-68.5 

-78.2 

3C 345 

4.00 

3.67 

0.92 

0.421 

0.406 

0.96 

- 

- 

- 

- 


0.702 

0.648 

0.92 

2.00 

2.07 

1.0 

5.72 

5.66 

-86.2 

-86.4 

NRAO 512 

0.887 

0.818 

0.92 

0.537 

0.628 

1.2 

- 

- 

- 

- 

SgrA* 

1.01 

0.946 

0.94 

2.58 

2.56 

0.99 

- 

- 

- 

- 


Column designation; 1 - source name; 2-3 - model flux density of the component[Jy]; 4 - ratio of the model flux density of Daejeon 
correlator to that of the DiFX; 5-6 - size of the component [mas]; 7 - ratio of the model size of the Daejeon correlator to that of DiFX; 
8-9 - angular distance from the central component [mas]; 10-11 - position angle[°]; 


observations in various observing modes (e.g., full track 
imaging mode). Since the comparison in this paper 
used the observations with short integration times of 
10-20 min (75min for SgrA*), we have an uncertainty 
in imaging results. The image quality values of of 
the final images are in the range of 0.53-0.88. 

4.2. Fringe Tracking 

We found that part (5%) of the total difference of 
8% in flux density between the two correlators are 
caused by other reaso ns which are not rel ated with the 
double-layer pattern. Ilguchi et. all (|2000ll investigated 
the performance of the hardware correlator developed 
for the VLBI Space Observatory Programme (VSOP) 
and found that the loss of visibility amplitude due to 
the hardware phase tracking is 4%. The loss almost 
corresponds to the 5% difference in flux density of im¬ 
ages using the double-layer-pattern free data (i.e., data 
for IF 8 and 9). For the Daejeon correlator, the fringe 


tracking is done in VCS in the same way as described in 
Ilguchi et. all(l2000ll . So a major part of the 5% degrada¬ 
tion in hardware correlation may come from the differ¬ 
ence in the way of fringe phase tracking. The additional 
1% of the difference can be attributed to other reasons 
(e.g., bit-jump correction). 

5. Conclusions 

We evaluated the performance of a new hardware cor¬ 
relator in Korea, the Daejeon correlator, by conducting 
KVN VLBI observations at 22 GHz. We correlated the 
acquired data with DiFX and the Daejeon correlator for 
comparing the visibility data and imaging results. The 
flux densities and brightness distributions of the tar¬ 
get sources for the two correlators are consistent with 
each other within < 8%. The difference is comparable 
with the amplitude calibration uncertainties of KVN 
observations at 22 GHz. We also found that the 8% 
difference in flux density is caused mainly by (a) the 
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difference in the way of fringe phase tracking between 
the DiFX software correlator and the Daejeon hardware 
correlator, and (b) an unusual pattern (a double-layer 
pattern) in the amplitude correlation output from the 
Daejeon correlator. The visibility amplitude loss by the 
double-layer pattern is as small as 3%. We conclude 
that the new hardware correlator produces reasonable 
correlation outputs for continuum observations which 
are consistent with the outputs from the software cor¬ 
relator, DiFX. 

Acknowledgments 

We thank Walter Alef, David Graham, and Alan Whit¬ 
ney for their fruitful discussion on the performance of 
the Daejeon correlator. We would like to thank the 
anonymous referee for important comments and sugges¬ 
tions which have enormously improved the manuscript. 
We are grateful to all staff members in KVN who helped 
to operate the array and to correlate the data. The 
KVN is a facility operated by the Korea Astronomy and 
Space Science Institute. The KVN operations are sup¬ 
ported by KREONET (Korea Research Environment 
Open NET work) which is managed and operated by 
KISTI (Korea Institute of Science and Technology In¬ 
formation) . 

Reeerences 

Bare, C., Clark, B. G., and Kellermann, K. I. 1967, 
Interferometer Experiment with Independent Local 
Oscillators, Science, 157, 189 
Carlson, B. R., Dewdney, P. E., Burgess, T. A., et 
al, 1999, The S2 VLBI Correlator: A Correlator for 
Space VLBI and Geodetic Signal Processing, PASP, 
111, 1025 

Casse, J. L. 1999, The European VLBI Network MkIV 
Data Processor, New A Rev., 43, 503 
Clark, B. G., Cohnen, M. H., and Jauncey, D. L. 1967, 
Angular Size of 3C 273B, ApJ, 149, L151 
Clark, B. G., 1973, The NRAO tape-recorder interfer¬ 
ometer system, Proc. IEEE, 61, 1242 
Deller, A. T., Tingay, S. J., Bailes, M., et al. 2017, 
DiFX: A Software Correlator for Very Long Baseline 
Interferometry Using Multiprocessor Computing En¬ 
vironments, PASP, 119, 318 
Deller, A. T., Brisken, W. F., Phillips, C. J., et al. 
2011, DiFX-2: A More Flexible, Efficient, Robust, 
and Powerful Software Correlator, PASP, 123, 275 
Greisen, E. W. 2009,The FITS Interferometry Data In¬ 
terchange Convention, AIPS Memo, 114 
Horiuchi, S., Kameno, S., Nan, R., et al. 2000, Imag¬ 
ing Capability of the Mitaka VSOP Correlator, Adv. 
Space Res., 26, 625 

Iguchi, S., Kawaguchi, N., Murata, Y., et al. 2000, De¬ 
velopment and Performance of the Real-Time VLBI 
Correlator (RVC), lEICE Trans. Commun., E83-B, 
2527 


Kondo, T., Koyama, Y., Nakajima, J., et al. 2003, In¬ 
ternet VLBI System Based on the PC-VSSP (IP- 
VLBI) Board, in ASP Conf. Ser. 306, New Technolo¬ 
gies in VLBI, ed. Y.C. Minh (San Francisco: ASP), 
205 

Lee, S.-S., Byun, D.-Y., Oh, C. S., et al. 2011, Single- 
Dish Performance of KVN 21 m Radio Telescopes: 
Simultaneous Observations at 22 and 43 GHz, PASP, 
123, 1398 

Lee, S.-S., Petrov, L., Byun, D.-Y., et al. 2014, Early 
Science with the Korean VLBI Network: Evaluation 
of System Performance, AJ, 147, 77 
Lobanov, A. P., Krichbaum, T. P., Witzel, A., et al. 
2006, Dual-Frequency VSOP Imaging of the Jet in 
S5 0836-f710, PASJ, 58, 253 
Moran, J. M., Crowther, P. P., Burke, B. F., et al. 1967, 
Spectral Line Interferometry with Independent Time 
Standards at Stations Separated by 845 Kilometers, 
Science, 157, 676 

Napier, P. J., Bargri, D. S., Clark, B. G., et al. 1994, 
The Very Long Baseline Array, Proc. IEEE, 82, 658 
Petrov, L., Lee, S. S., Kim, J., et al. 2012, Early Science 
with the Korean VLBI Network: The QCAL-1 43 
GHz Calibrator Survey, AJ, 144, 150 
Rogers, A. E. E., Cappallo, R. J., Hinteregger, H. F., 
et al. 1983, Very-long-baseline radio interferometry 
- The Mark HI system for geodesy, astrometry, and 
aperture synthesis, Science, 219, 51 
Ryle, M., and Hewish, A. 1960, The synthesis of large 
radio telescopes, MNRAS, 120, 220 
Shepherd, M. C., Pearson, T. J., and Taylor, G. B. 
1994, DIFMAP: an interactive program for synthesis 
imaging, BAAS, 26, 987 

Shibata, K. M., Kameno, S., Inoue, M., et al. 1998, Mi¬ 
taka Correlator for the Space VLBI, in ASP Conf. 
Ser. 144, Radio Emission from Galactic and Extra- 
galactic Compact Sources, ed J.A. Zensus et al. (San 
Francisco: ASP), 413 

Thompson, A. R. 1999, in Synthesis Imaging in Radio 
Astronomy II, eds. G. B. Taylor, C. L. Carilli, & R. 
A. Perley, ASP Conf. Ser., 180, 11 
Whitney, A. R. 1993, The Mark IV VLBI Data- 
Acquisition and Correlation System, in lAU Symp. 
156, Developments in Astrometry and Their Impact 
on Astrophysics and Geodynamics, ed. 1. 1. Mueller 
and B. Kolaczek (Dordrecht:Kluwer), 151 
Whitney, A. R. 2002, Mark 5 Disc-Based Gbps VLBI 
Data System, in Proc. 6th European VLBI Net¬ 
work Symp., New Developments in VLBI Science and 
Technology, ed.E. Ros et al. (Bonn: Max-Planck- 
Institut fiir Radioastronomie, 41 
Wietfeldt, R., Baer, D., Cannon, W. H. et al. 1996, The 
S2 very long baseline interferometry tape recorder, 
IEEE Trans. Instrumentation and Measurement, 45, 
923 

Wilson, W., Roberts, P., and Davis, E. 1996, in 
Proc. 4th ATP workshop, ed. E. A. King (Sydney: 
ATNF/CSIRO), 16 



